BNL-HET/06/6 
KANAZAWA-06-09 
KEK-CP-177 



Vector form factor in Ki 3 semileptonic decay with two flavors of 
dynamical domain-wall quarks 

Chris Dawson 1 , Taku Izubuchi 1,2 , Takashi Kaneko 3,4 Shoichi Sasaki 1 ' 5 , and Amarjit Soni 6 
1 RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 
2 Institute for Theoretical Physics, Kanazawa University, Kanazawa, Ishikawa 920-1192, Japan 
A High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 
4 Graduate University for Advanced Studies, Tsukuba, Ibaraki 305-0801, Japan 
J Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 
6 Physics Department, Brookhaven National Laboratory, Upton, NY, 11973, USA 

(Dated: February 7, 2008) 

We calculate the vector form factor in K — > ttIu semileptonic decays at zero momentum transfer 
f+(0) from numerical simulations of two-flavor QCD on the lattice. Our simulations are carried out 
on 16 3 x 32 at a lattice spacing of a ~ 0.12 fm using a combination of the DB W2 gauge and the 
domain- wall quark actions, which possesses excellent chiral symmetry even at finite lattice spacings. 
The size of fifth dimension is set to L s = 12, which leads to a residual quark mass of a few MeV. 
Through a set of double ratios of correlation functions, the form factor calculated on the lattice 
is accurately interpolated to zero momentum transfer, and then is extrapolated to the physical 
quark mass. We obtain /+(0) =0.968(9) (6), where the first error is statistical and the second is the 
systematic error due to the chiral extrapolation. Previous estimates based on a phenomenological 
model and chiral perturbation theory are consistent with our result. Combining with an average 
of the decay rate from recent experiments, our estimate of /+(0) leads to the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element |Vu S | =0.2245(27), which is consistent with CKM unitarity. These 
estimates of /+(0) and \V US \ are subject to systematic uncertainties due to the finite lattice spacing 
and quenching of strange quarks, though nice consistency in /+ (0) with previous lattice calculations 
suggests that these errors are not large. 



I. INTRODUCTION 

There has recently been rapid progress in the precise determination of the elements of the CKM 
matrix [l| leading towards a stringent test of its unitarity. Let us recall that such a test is a 
powerful method to search for new physics beyond the standard model. In particular d — > u and 
s — > u semileptonic transitions provide the most precise constraints on the size of the elements, 
and hence CKM unitarity on the first row 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1-5 (1) 

can now be examined accurately @, [1[ . 
The values quoted in the 2004 PDG Q 

|K<j| = 0.9738(5), | K s | = 0.2200(26), \V ub \ = (3.67 ± 0.47) x 10~ 3 , (2) 



lead to 

5 = 0.0033(15), (3) 

which deviates from zero by two a. We have to improve the accuracy on 5 in order to confirm 
whether this deviation is a genuine signal of unitarity violation. We note that \V u b\ is so small that 
it can be safely neglected in this unitarity test. Since about half the error of (5 comes from the 
uncertainty in |V^d|, and another half from |V^ S |, we need a more precise determination of both of 
these two elements. In the present paper, we focus on the determination of |V^ S |. 
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We also note that \V US \ gives the basic parameter A in the Wolfenstein parameterization of the 
CKM matrix A precise determination of \ V US \ is, therefore, important also for phenomenological 
studies of CP violation processes based on this parametrization. 

So far, \V US \ has been determined from several processes: Ki 3 decays 0| which provide the value 
in Eq. ©, hyperon j3 decays @, and 7T M 2 decaysQ, and hadronic r decaysQ. At the moment, 
the K13 decays provide the most precise determination among these, and its result has been quoted 
in the PDG. We, therefore, try to determine \V US \ through the K13 decays. 

As will be explained in SecHH \ V US | can be determined through experimental determination of the 
decay rate T and theoretical calculation of the vector form factor at zero momentum transfer .f+ (0) ■ 
The two a deviation from unitarity in Eq. ([3]) motivated recent measurements of T [9I. [Tol [TT1. fl3j . 
These experiments prefer a larger value for \V US \, which is consistent with unitarity. However, a 
precise calculation of /+(0), say with an accuracy of 1 %, is also indispensable in order to establish 
this consistency with unitarity. 

A good theoretical control on /+(0) is provided by SU(3) symmetry and the Ademollo-Gatto 
theorem (l4j . The vector current conservation guarantees /+(0) = 1 at zero momentum transfer 
[lH . and then the Ademollo-Gatto theorem states that the SU(3) breaking effects /+(0)— 1 start at 
second order in (m s —m u d), where m u d represents the averaged mass of up and down quarks. This 
also indicates that the leading correction to /+(0) in chiral perturbation theory (ChPT) does not 
contain the low-energy constants (LECs) of the next-leading order chiral Lagrangian, and hence is 
practically free of uncertainties. 

In this paper, we calculate the vector form factor of the K13 decay at zero momentum transfer 
/+(0) non-perturbatively from numerical simulations of lattice QCD with two degenerate flavors 
of dynamical quarks, which are identified with up and down quarks. Strange quarks are treated in 
the quenched approximation. To make the best of use of the good theoretical control mentioned in 
the previous paragraph, we employ a combination of the DBW2 gauge action and the domain- 
wall quark action [13, Oil' which has excellent chiral symmetry even at finite lattice spacings. We 
also employ the so-called double ratio method [n| to improve the accuracy on the form factor. 
Preliminary results of these calculations have been reported in Ref . [20] . 

This paper is organized as follows. We present a brief introduction of the K13 decays and status 
of experimental and theoretical studies on them in Sec. ITT] Our simulation method is introduced in 
Sec. IIIII Section IPV1 is devoted to our determination of the form factor at finite momentum transfer 
from a double ratio of three-point functions. We describe the interpolation of the form factor to 
zero momentum transfer and the chiral extrapolation in Sees. fVl and IVII Section IVlII presents our 
estimate of \V US \. Our conclusions are given in Sec. IVIIll 



II. K l3 DECAYS 
A. phenomenology of Kn decays 

The K13 decays are K to n semileptonic decay channels 

K° 3 : K° -> tt-I+h, (4) 
K+ : K + ^-k°1 + v u (5) 

where I represents the electron or muon. In the following, we mainly consider the neutral kaon 
decay K® s . A simplification in theoretical studies of these decays is that the matrix element of 
the axial current vanishes due to the parity symmetry. Therefore, the decay amplitude contains 
only the matrix element of the vector current — sj^u which can be expressed in terms of form 
factors 

(n(p')\V,\K(p)) = (p,+p',) f + {q 2 ) + (p»-p',) f-(l 2 ), (6) 
where q~p — p' represents the momentum transfer. 



3 



In literature, the so-called scalar form factor 



f Q (q 2 ) = f + (q 2 ) + M /_ M J -(<1 2 ), (7) 



and 



are often used instead of f-{q 2 ). In particular, fo(q 2 ) is a useful quantity in lattice calculations, 
since i) it equals to /+(0), which appears in the expression of the decay rate (see Eq. (fTO)) below), 
at zero momentum transfer q 2 = 0, ii) it can be precisely calculated from the matrix element with 
kaon and pion momenta equal to zero 

(7r(0)|^|ir(0)) = (M K +M n )f (q 2 max ), (9) 

where q^ ax = (M K -M W ) 2 . 

The rate of the Ki 3 decays is given by [15| 

G 2 

r = t^m 5 k c 2 i\v us \ 2 \u(o)\ 2 s c M + s c1 u), (10) 

where / and S cw (l +5 em ) represent the phase space integral and radiative corrections, respectively. 
The Clebsch-Gordan coefficient C 2 = 1(1/2) for the neutral (charged) kaon decay is written ex- 
plicitly in the above expression so that /+(0) for both decay channels equals to unity in the SU(3) 
symmetric limit. 

The phase space integral I is generally defined as 

U(q 2 ) 3M 2 (M 2 K -M 2 ) 2 Mq 2 ) \ 

f + (q 2 ) + {2q 2 + M 2 )\ I U(q 2 )j' ( ' 

Xi = q 4 + M$ ( + M*-2q 2 M 2 K -2q 2 M%-2M%M%, (12) 

where qo is a reference value of the momentum transfer. If we take go 7^ 0, /+(0) in Eq. (fT0|) has 
to be replaced by f+(qo). As in Eq. (fTU|) . q is usually set to so that / depends on f+,o{q 2 ) only 
through small coefficients A+\ A+ and which parametrize the q 2 dependence of / +i o(? 2 ) 

f+(q 2 ) = /+(0) (l + A^^ + Afg 4 ), (13) 
fo(q 2 ) = /o(0) (l + A^g 2 ), (14) 

where we include the quadratic term suggested by the KTeV experiment to f+(q 2 ) [2l|. From 
recent experimental measurements of these coefficients [111, [22|, [23[ , the current estimate for / is 
- 0.154 (0.159) for K° 3 (K+ 3 ) and 0.102 for K° 3 with an accuracy of around 1% @,S], where a 
dominant error comes from the choice of the parametrization form of the q 2 dependence of f+(q 2 ). 
The choice of the reference scale qo=0 forces us to study the q 2 dependence of the form factor and 
take the limit of q 2 =0 as in Sec. |V] 

The radiative corrections split into the short-distance electroweak piece S cw and the long-distance 
electromagnetic piece (1 + S em ). The former is precisely determined as S cw — 1.022 [24j |. Chiral 
perturbation theory including the electromagnetic interaction [25| and a phenomenological model 
|2a | reveal that the latter is small correction and its uncertainty leads to < 1% error to T. 

Recently several new experimental determination of T have been performed to clarify the origin 
of the two a deviation from unitarity in Eq.© 0, [TFj, [TJ, [TJ, [TJ]. In Ref.@], \V US /+(0)| = 0.2173 
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is obtained from the new measurements of T and recent estimates of I and <5 om through Eq. (| 10|) . 
This is about 3% larger than 0.2114 corresponding to \V US \ in Eq. J2]), and may lead to a good 
consistency with CKM unitrarity. However, we have to determine /+(0) with an accuracy of 1% 
in order to make a definite conclusion on this unitarity test. 



B. previous theoretical studies of /+(0) 

In previous theoretical studies, /+(0) is considered in the following ChPT expansion 

/+(0) = 1 + fa + h + 0(p 6 ), (15) 

where fa n is the 0(M 2n K ) — 0(m,g) correction to /+(0). We note that the leading term is unity, 
because /+(0) becomes the Clebsch-Gordan coefficient in the SU(3) symmetric limit thanks to the 
vector current conservation, and we explicitly factor it out Eq. ()10l) . 

The Ademollo-Gatto theorem (T^| states that SU (3) breaking effects are second order in (m s — 
Tn u d)- From previous theoretical studies, it turned out that the SU(3) breaking effects /+(0) — 1 
are order of 3-5%. Therefore, we can achieve 1% accuracy on /+(0) by calculating the SU(3) 
breaking effects with an accuracy of 20-30%, which is not prohibitively challenging. 

In addition, the Ademollo-Gatto theorem guarantees that the leading correction fa does not 
contain any poorly-known LECs, which are associated with analytic terms from the 0(p 4 ) chiral 
Lagrangian. For example, its ChPT formula for the Kf 3 decay is given by [27j 

1 3 

fa = h k»tv + 2 H k+it + 2 H k+v + tV3{H K iT - H Kv ), (16) 



i r „ „ iMiMi 



Urn = { .\fr, + .Ur,+ r In 



IPQ = _ 64tf2f% I ^ ~Mp — Ml 



M 2 



(17) 



where e= (\/3/4) (m d - m»)/(m, - m ud ). 

However, the next leading correction fa contains LECs in 0(p 4 ) and 0(p 6 ) chiral Lagrangians 
[H, US H3 [U, H3- Therefore, fi is difficult to determine only from ChPT, and hence a phe- 
nomenological estimate /t = — 0.016(8) by Leutwyler and Roos [15| has been employed in previous 
determinations of \V US \. 

Clearly, it is desirable to calculate /+(0) non-perturbatively. This background led to the first 
lattice study in quenched QCD [33]. They demonstrated that lattice calculations can achieve the 
1% accuracy for /+(0) by using a set of the so-called double ratios of correlation functions, and 
by making good use of the ChPT formula for /2, namely Eq. (|16p . in the chiral extrapolation of 
their lattice data. They employed the non-perturbatively 0(a)-improved Wilson quark action and 
obtained /+(0) = 0. 960(9), which is consistent with the Leutwyler-Roos's estimate. 

The calculation was extended to two-flavor QCD with the (9(a)-improved Wilson quark action 
[34| and to three-flavor QCD with an improved Kogut-Susskind (KS) quark action [35|]. While the 
q 2 dependence of f+(q 2 ) has not been investigated in the latter study, their estimates of /+(0) are 
consistent with that in quenched QCD. 

It is also worth while to note that the so-called twisted boundary condition (36|, which enables 
us to explore small q 2 region, has been tested in quenched QCD [37J. 



III. SIMULATION METHOD 



In this study, we calculate the kaon form factor /+(0) by numerical simulations of lattice QCD 
with two-flavors of dynamical quarks, which are identified with up and down quarks. We note 
that the isospin breaking effects in Eq. (fTo]) is proportional to e w 0.01, and fa itself is of order 
2-3% shift, so the correction due to isospin breaking is well below our target accuracy on /+(0) 
of 1%. Strange quarks are treated in the quenched approximation. We employ the domain-wall 
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quark action [17| , which has the following advantages over the conventional Wilson- and KS-type 
actions. 

First, it possesses chiral symmetry even at finite lattice spacings in the limit of L s — > oo, while 
the conventional Wilson- and KS-type fermions break the chiral symmetry explicitly. The chiral 
behavior of physical quantities may be distorted with the conventional fermions. Rigorously speak- 
ing, we have to take account of effects of the explicit symmetry breaking in the chiral extrapolation 
of the quantities obtained with the conventional fermions at finite lattice spacings [38l . [39| . In con- 
trast, chiral extrapolations with domain- wall quarks are fairly straightforward and simple as they 
need use essentially continuum ChPT [40| • This is particularly important in this study of the K13 
from factor, since we can safely subtract the leading correction f 2 from our lattice data of /+(0) by 
using the ChPT formula Eq. lfTB)) before the chiral extrapolation so that systematic uncertainties 
due to the extrapolation influence the final result only through the small higher order corrections. 

Another advantage of the use of domain-wall quark action is that it is automatically 0(a)- 
improved. Unlike the Wilson-type fermions, the non-perturbative tuning of improvement coeffi- 
cients for the vector current is not necessary to remove possibly large O(a) effects [4l|. We also 
note that the leading 0(a 2 ) scaling violation in physical quantities such as the kaon B parameter 
is not large at a" 1 ~ 2 GeV @,|43j]. 

These virtues of domain-wall fermions are lost if L s is not sufficiently large. Since the CPU cost 
is proportional to L s , it is limited to values around 10-20 in practical unquenched simulations. In 
order to improve chiral properties of the form factors with L s fixed, we employ the DBW2 gauge 
action [l6[ with which the light hadron spectrum and the kaon B parameter show better chiral 
properties than with the conventional plaquette action [44l |45j. 

We use gauge ensembles generated on a i 3 x T — 16 3 x 32 lattice at f3 = 0.80, as discussed in 
Ref . [45| . The lattice spacing is a ~ 0.12 fm and the physical spatial size is La ~ 1.9 fm. We set the 
domain-wall height to M$ = 1.8 and the fifth-dimensional length to L s = 12. The resulting residual 
quark mass is a few MeV. We simulate sea quark masses m U( i S oa = 0.02, 0.03, and 0.04 in the range 
of m S! phy S /2 < w«d.sca ^ ^s,phys, where m S!P h ys represents the physical strange quark mass. Our 
statistics are 94 configurations separated by 50 HMC trajectories at each sea quark mass. We refer 
to Ref. [!|| for further details on the configuration generation. 

On these gauge ensembles, we calculate two and three point functions 



C p (t;p) = ^2(o PtSnk (x,t + to)0^ src (0,t ) 



> Z hrc Z P,snk Ep(p)t 

t^oo 2E P (p) 1 ' 

C^(t,t';p,p') = £ (o Q , snk (x' ) t' + 1 ) V„(x, t + t )O PiSIC (0, to)) e^*'-*)^, 



X, X' 



_Zp^ic_^Qfn^__ 

*,(t'^)-oc ; 4E P {p)E Q (p')Z v ^" 



(Q(p')\vW\P(p)) 

xe -Er(p)t-E Q (p')(t'-t) ^ 



where P and Q denote K or 7r meson. The sink (source) operator for the meson P is represented 
b y C 'p, ) snk(src)' and its over l a P t0 tnc physical meson state is given by Z P ^ snk{src) = (0|Op !Snk(src) |F). 
We denote the energy of meson P with a spatial momentum p by Ep(p). The renormalized vector 
current with the renormalization factor Zy is represented by . 
In our preliminary study [201 ] , an exponential smeared operator 

£>(|r|)g(x) 75 <7(x + r), cj>{\v\) = A exp [-B |r|] (20) 



with A=l.2 and B = QA was used for the initial meson. We observe that the correlators C p (t; p) 
and C p ®(t, t'\ p, p') with non-zero p show poor signals with this choice of the smearing function. 
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This may suggest that this operator is too close to the wall source, and the correlators have 
small overlap with meson states with non-zero p. In this study, therefore, we use more localized 
operators with B = 0.5, 0.6, and 0.7 at m u d = 0.02, 0.03, and 0.04, respectively. We take a single 
choice of ^o = 4 for both of C p (t;p) C P ®(t, t'\ p, p'). The local operator, which is combined with 
the sequential source method [46[ for C p ®(t,t';p, p'), is used for the sink meson operator. 

As in Ref.[45[, we calculate the quark propagator with each of the periodic and anti-periodic 
boundary conditions in temporal direction for quarks. The correlation functions C p and C P ® are 
constructed by using the averaged quark propagator over the boundary conditions. This procedure 
cancels effects of valence quarks wrapping the lattice in the temporal direction by odd number of 
times, and enable us to take the time slice for the sink operator (t'+to) for C P ® larger than T/2. 
In this study, we fix t' to 24 for all combinations of sea and valence quark masses. 

In our measurement of C p and C P ® , we fix the valence ud quark mass equal to the sea quark 
mass, and take four strange quark masses m s =0.02, 0.03, 0.04, and 0.05, which are roughly in a 
range of [m SiP h ys /2, 5 m SiP h ys /4]. For the meson momentum, we take all possible configurations p 
with |p| 2 = 0, 1, and 2 for the initial meson "P" in C p and G P ® . For notational simplicity, we use 
the momentum p% in units of 2-7T / 'Lk throughout this paper. Two configurations p = (0,0,0) and 
(—1,0,0) are used for the final meson "Q" in C P ® . 



IV. SCALAR FORM FACTOR AT qi 



A. Meson masses 



In extraction of the form factors from C P ® , we need precise knowledge of the pion and kaon 
masses and their energies with finite momenta, which appear in Eq. (|19p . In Fig. [TJ we plot 
the effective mass for pions calculated from C n . We observe a clear and long plateau in the 
effective mass, and hence the (lattice) meson masses summarized in Table U are determined with 
an accuracy of < 1%. We note that masses presented in Ref . (45| . in which the local and wall 
sources are employed, are consistent with ours within statistical error. 

For the meson energies with non-zero momenta, we use an estimation from the fitted mass M 
and the lattice dispersion relation 



£(P) 2 



M 



E(p) = 2sinh 



E(p) 



Pk 



2 sin 



Pk_ 
2 



(21) 



instead of the fitted energy to C p , since its statistical error rapidly increases as the size of the 
meson momentum increases. We observe that i) the fitted energy shows a good consistency with 
Eq. (|2"T|) as in Fig. O ii) the final result for /+(0) at the physical quark mass does not change 
significantly if we use the dispersion relation in the continuum limit instead of Eq. (|21[) ■ 



B. /o(<?max) 

We consider the double ratio which was originally proposed in Ref. [19] for B meson decays 

m .s _ cf*(t,^o,o)Q*(M';Q,o) (99 , 

H[t) Cf K (t,t';0,0)Cr(t,t';0,Oy 

where we fix t' to 24 as mentioned in Sec. lIII[ and t' dependence of R is ignored in the following. All 
of Zy, Zp iSIC , ^Q jSn k and the exponential damping factor in Eq. (|19[) are exactly canceled in R(t). 
As a result, R(t) contains only meson matrix elements with zero momentum, and gives /o(<7max) 

R{t) (K\vt ] \K){.\vt ] \-) ~ * M * M « ' 
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We note that R(t) is exactly equal to unity in the SU(3) symmetric limit and hence it is a useful 
quantity to measure SU(3) breaking effects to /o(<?max)- 

In Fig. [31 we show the three-point functions C^* , C\ K , Cf K and C%* for each jackknife sample. 
Their fluctuation leads to the jackknife error of about 5%. As seen in the figure, they are highly 
correlated with each other, and hence the double ratio R(t) has a very small fluctuation over the 
jackknife samples, which leads to the error of about 0.03%. 

Figure [5] shows that the magnitude of the statistical error does not change drastically at other 
values of t. From the fitted value of R(t) and measured meson masses on lattice, we obtain /o(<?max) 
summarized in TablcHIl There is a tendency that the error of /o(<Zmax) increases as m s deviates from 
m u d, which is probably because /o(<7max) deviates from its trivial value 1 towards larger \m s —m u d\- 
However, the statistical accuracy is < 0.1% even in the worst case (m u d, m s) = (0.02, 0.05). 



V. INTERPOLATION TO ZERO MOMENTUM TRANSFER 

To study the q 2 dependence of the form factor, we calculate 

F(p,p) = y-r^ r 1 + . , , » C(g ) , (24) 

JO «ax ) V E K (p) + K (p'j / 



from a ratio 

i?(t;p,p') 



C^(t,ti;p,pi)C K (t;0)C"(t> -t;0) 
Cf »(t, t>; 0, 0) C*(t; p) C*(tf - i; p') 
£?jf(p)+^(p' 



f,(t'-t)— >oo Mr- + Mjr 



F(p,p'). (25) 



Since R(t; p, p') has its trivial value 1 at |p| = |p'| = 0, it might be a good probe to study how 
fo(q 2 ) changes as q 2 deviates from g^ ax . We also note that R(t;p, p') is reduced to the double 
ratio employed in Ref. [HJ|, if we fix p = 0. 

In order to reduce the statistical error, we calculate R(t;p,p') from three-point functions aver- 
aged over momentum configurations which correspond to the same momentum sizes |p| and |p'|. 
Figure [5] shows example of F(p, p') obtained from R(t; p,p') as a function of t. We observe that 
data with the smallest non-zero momentum |p| 2 = I show a clear plateau, and hence F(p,p') is 
determined with an accuracy of roughly 5%. 

In order to convert F(p, p') to fo(q 2 ), we evaluate £,(q 2 ) by employing the method proposed in 
Ref. ,33|. Namely, we measure the double ratio 



R (f n W) C^(t,t>;p,pi)C« K (t,t>;p,p>) 



and calculate ) from 



p{ 2, = -(E K (p) + E K (p>)) (p+p% + (E K (p)+E 7T (p 1 )) (p+p% Rk . s 

m> {E K {p)+E K {pi)){p-p>) k -{E K {p)-E^p>)){p+p>) k R k ■ {l) 

As for R(t; p, p'), we first take the average for the relevant three-point functions over appropriately 
chosen momentum configuration and the Lorentz index for (/x=l,2,3), and then double ratio 
Rk(t;p, p') is constructed from the averaged correlation function. We note that Rk(t;p, p') is 
exactly unity in the SU(3) symmetric limit, and is sensitive to SU(3) breaking effects in the 
matrix element . 

Figure [5] shows examples of R k (t;p, p') as a function of t. We observe that, at most of our 
simulated quark masses, R k (t;p,p') is close to unity, and hence £,(q 2 ) from Eq. (f2"T|) has small 
magnitude <0.1. Its error is typically 30% -100% with our statistics. 
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We note that F(p, p') and £(g 2 ) can be calculated also from the following ratios constructed 
from the tt — > K matrix element 

5 V , „ = CZ K (t,t';p,p')C*(t;Q)C K (ti-t;0) 
l ' P ' PJ CZ K (t,t';0,0)C*(t;p)C K (t>-t;p'y { ' 



; • 1 1 ' " C][«(t,t>;p,p>)C™(t,t>;p,p>y U1 



We confirm that, for |p| = and |p'| = 1, R(t;p,p') and R'(t;p',p) give consistent results for 
F(p, p'), while the latter leads to much smaller error. This is because, as described in Sec JIIIl 
C^(t,t';p,p') is measured with the single choice of the final meson momentum p' for each |p'|, 
and we can not take average of R(t;p,p') over the momentum configuration "{p'}". We also 
observe that data with |p|, |p'| >0 show poor signal. Therefore, in the following analysis, we use 
F(p, p') and £(q 2 ) obtained from R(t; p, p'), Rk(t; p, p'), R'(t; p, p'), and R' k (t; p, p') with |p'| =0. 

In order to take the limit of zero momentum transfer of fo(q 2 ) reliably, we test two methods to 
calculate /o(0), and check the consistency between results from these methods. In the first method, 
which we call method-1 in the following, we calculate fo(q 2 ) at simulated q 2 from F(p,p') and 
£,{q 2 ). Then, the results for fo(q 2 ) and /o(<Zm ax ) horn R(t) are interpolated to q 2 = 0. As well as 
the linear fitting form Eq. (fT4|) , we test the quadratic form 



and the pole form 



fo(q 2 ) = MOMl + A^V + A<V), (30) 



fo(q 2 ) = / (0)/(1-A<V). (31) 



The physical value of qf nax for decays is much smaller than that for B meson decays. In 
addition, (?^ ax is further reduced in our lattice calculation, since the simulated values of m u d 
are larger than its physical value. As shown in Fig. /o(0) can be determined by a very short 
interpolation from <7 2 lax , where we have very accurate data /o(<Zmax) from R(t). As a result, the 
choice of the interpolation form does not affect the interpolated value /o (0) significantly. Actually 
we observe that the interpolated values obtained from the three forms, Eqs. (fH)) . (|50"|) and (f5T|) 
are consistent with each other. 

In the following, we employ the result from the pole form Eq (|3Tj) . because i) data of fo(q 2 ) does 

(2) 

not have a strong curvature, and hence the quadratic fit leads to ill-determined Aq , which has 
typically 100% error, and ii) the pole form leads to slightly smaller value of x 2 /dof than the linear 
fit. The fit parameters are summarized in Table HTT1 While the pole form at m u d — 0.03 leads to 
a slightly higher value of x 2 /dof than at other values of m u d, the inclusion of the quadratic term 
does not reduce x 2 /dof significantly. We note that the statistical error on /o(0) is < 0.3% level. 

We also test an alternative method to calculate /o(0) employed in Ref. [34[ , which we call method- 
2 in the following. In this method, we first take the limit of F(p, p') and £(g 2 ) to q 2 =0, and then 
calculate /o(0) from F(p, p')|^2 =0 and £(0). Since i^(p,p') depends on two momenta |p|, and |p'|, 
the q 2 interpolation of F(p,p') has to be carried out using data with fixed |p| (or |p'|). This also 
enables us to identify |p'| (|p|) corresponding to <z 2 =0, which is needed to convert F(p,p')| g 2 =0 to 
/o(0). In the following, we repeat the interpolation for two data sets with |p| =0 and |p'| =0, and 
take the average of results for F(p, p')l g 2 =0 - 

For the q 2 interpolation, we test linear, pole, and quadratic fitting forms similar to Eqs. HH), 
and (|3"TT) , and employ the quadratic fit 

F(p,p% 2 = F(p,p')\ g2=0 -(l + c F:1 q 2 + c Ft2 q 4 ) (|p| or |p|' is fixed), (32) 



since this fit leads to the smallest value for % 2 /dof among the tested forms, and to the reasonably 
well-determined cf.i- Examples of this quadratic fit are shown in Fig. [5] 
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While £(0) has to be determined by an extrapolation, we observe that £,{q 2 ) has very mild q 2 
dependence, as seen in Fig. [9l and the simplest linear form 

H(q 2 ) = £(0)(l + c C ,i<7 2 ) (33) 

leads to reasonably small x 2 /dof. Table HVl shows £(0) and /o(0) calculated from F{p,p')\ q2=0 
and £(0). 

We observe that two methods of the q 2 interpolation give consistent results for /o(0) with each 
other. However, the error from method-2 is slightly larger than from method-1 due to the error for 
£(0) enhanced by the extrapolation. Therefore, we employ results from method-1 in the following. 

We note that the above observation is opposite from that in our preliminary study [20| , where 
method-1 give larger error for /o(0) mainly due to the uncertainty in £(q 2 ) at simulated q 2 . As 
described in Scc lIIIl we calculate correlation functions with the different choice of the smearing 
function from that in Ref.[20[. This improves accuracy in £(g 2 ) at simulated q 2 and hence /o(0) 
in method-1. However, in order to reduce the uncertainty from method-2, the change of the 
smearing functions is not sufficient and we need to have data at small q 2 for a better control of 
the extrapolation of £(q 2 ). 



VI. CHIRAL EXTRAPOLATION 



In the chiral extrapolation of /+(0) = /o(0), we rewrite the ChPT expansion Eq. fT5]) as 

/+(0) = l + / 2 +A/, (34) 

where A/ represents all higher order corrections starting at 0{M^ w ). 

As mentioned in Sec|TTl the leading correction fa does not have analytic terms from the 0{p A ) 
chiral Lagrangian thanks to the Ademollo-Gatto theorem. It is shown in Ref. [47| that the above 
statement is true even in two-flavor partially quenched (PQ) theory, as seen in their PQChPT 
formula of fa 

(PQ) = 2M 2 K + M 2 3M 2 K M 2 ln[M 2 /M 2 K ] 
12 32 7TV 2 64n 2 f 2 (M 2 K -M 2 ) 

M% (4 M 2 K - Ml) ln[2 - M 2 JM 2 K ] 

64tt2 f "2 {M\ -Ml) ■ 1 > 

Therefore, fa does not contain any poorly-known LECs in the 0{p A ) chiral Lagrangian, and its 
value at simulated quark mass can be precisely calculated from the measured meson masses through 
Eq. dSSJ). 

Consequently, the chiral extrapolation of /+(0) is nothing but the extrapolation of the higher 
order correction A/. Since it is also proportional to {m s — rn u( j) 2 thanks to the Ademollo-Gatto 
theorem, we consider the ratio 

RAf = {mi -Mir (36) 

as in Ref. 33], and extrapolate it to the physical quark mass. To this end, we test the following 
constant, linear and quadratic fits 

Raj = c , (37) 

RAf = c + c hv {M 2 K + M 2 ), (38) 

RAf = c + Cl . s M 2 + c hv {M 2 K + M 2 ), (39) 

RAf = c Q + c hs Ml + c hv {M 2 K + Ml)+c 2 , s M* + c 2 , v {M 2 K + M 2 ) 2 . (40) 
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The constant fit Eq. ([57]) should work if Af is dominated by the analytic term in f±. Linear and 
quadratic dependences in Eqs. (|38[) - (|40|) are assumptions for an effective description of the chiral 
logarithms in fi and higher order corrections. 

Figure [TU] shows the chiral extrapolation using the linear forms Eqs. ([55]) and (I39[) . We observe 
that i?A/ has mild dependence on the sea and valence quark masses, and the linear and even 
constant fits achieve sufficiently small value of x 2 /dof. While the quadratic fit Eq. (|40|) also gives 
a small value of x/dof, it leads to more than 100% error for both of C2, s and C2, v . We, therefore, 
do not use the results from the quadratic fit in the following discussion. 

From the fit parameters summarized in Table [V] and the physical meson masses determined in 
Ref. [45[ , we obtain Af at the physical quark mass which is also collected in Table [V] We note 
that all fits lead to consistent results for Af with each other. We obtain 

Af = -0.0Q9(9)(6), (41) 

by employing result from the linear fit Eq. (|39p . which is also employed in the unquenched calcu- 
lations in Refs.JU [Hj]. The first error is statistical, and the second is a systematic error due to 
the chiral extrapolation which is estimated as the largest deviation in Af among the constant and 
linear fits. 

By using /a = — 0.023 at physical quark mass in full QCD, we obtain 

/+(0) = 0.968(9)(6), (42) 

which is consistent with the previous lattice calculations [H, [34], [HI], employing different dis- 
cretizations (Wilson and KS) in quenched and unquenched QCD, listed in Table IVI[ as well as 
with estimates based on 0(p 6 ) ChPT [3(1 HH H2] , and the Leutwyler-Roos's value [l5l |. 

In Eq. (|42p , we have not included systematic uncertainties due to the discretization error and 
effects of dynamical strange quarks, which are difficult to estimate reliably without simulations 
at different lattice spacings or those in three-flavor QCD. However, these uncertainties affect our 
estimate of /+(0) only through the small higher order correction Af, and hence are expected not 
to be large. This is supported by the nice consistency with results from different lattice actions 
and/or with different numbers of flavors for dynamical quarks (see Table fVTl) . We also note that the 
RBC and UKQCD Collaborations have already started large-scale simulations with three flavors 
of dynamical domain- wall quarks, and their preliminary estimate is consistent with Eq. (|42[) . In 
particular, that comparison should provide a more reliable estimate of systematic error due to 
quenching of strange quarks. 

We also calculate £(0) at the physical quark mass. Since £(0) vanishes in the SU(3) symmetric 
limit, we test the following simple linear fit 

m=diAM 2 K -Ml), (43) 

and find that this leads to a reasonable value of x 2 /dof as shown in Table IVIII The fit line is 
plotted in Fig. QTJ We obtain 

£(0) = -0.105(22), (44) 
which is consistent with the experimental values —0.01(6) for K® 3 and —0.125(23) for Q]. 



VII. \V US \ AND CKM UNITARITY 

By combining with an estimate of |V^ S /+(0)| = 0.2173(8) based on the recent experimental 
determination of T |3j , we obtain 

|K,| = 0.2245(26) (8), (45) 

where the first and second errors come from the uncertainty in / + (0) and | V us f+ (0) | , respectively. 
This leads to 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1-5, 5 = 0.0013(16), (46) 
which is completely consistent with CKM unitarity. 
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VIII. CONCLUSION 



In this paper, we have calculated /+(0) from numerical simulations of two-flavor dynamical QCD 
using the domain-wall quark action. We obtained 

j+(0) = 0.968(9)(6), \V US \ = 0.2245(26)(8), (47) 

which supports CKM unitarity. While we have not estimated systematic uncertainties due to the 
use of the finite lattice spacing and the quenched approximation for strange quarks, these are 
expected to be small from the nice consistency with other lattice estimates. 

Our result for /+(0) is consistent with the phenomenological estimate, which has been used 
in previous determinations of \V US \, and hence has not changed \V US \ significantly. The main 
significance of this study is that now /+ (0) has been calculated non-perturbatively from two-flavor 
QCD and its uncertainties can be systematically reduced in future lattice calculations. Systematic 
errors which all the present lattice calculations share are those connected with the interpolation in 
momenta and extrapolation in mass. In both these cases it was necessary to use an ansatz. Since 
the interpolation in the momentum transfer was over a very small range, and the extrapolation in 
mass systematically took into account the calculated behavior up to NLO in ChPT, relying on the 
ansatz only for higher orders, both these effects are expected to be small. 

However, it should be noted that the sea quark masses used in this calculation are relatively 
heavy, and to be confident that ChPT is a good description of the data it would be advisable 
to move to smaller masses. In turn this will make the momentum interpolation more difficult 
as g^ ax deviates further from 0. Another important step in the future is clearly an extension to 
dynamical three-flavor QCD. The RBC and UKQCD Collaborations' study of three-flavor QCD is 
well underway [48|, |49(, and a more reliable estimate of systematic uncertainties in /+(0) will come 
in the near future. 
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TABLE I: Fitted meson masses. 



m u d 




M n 


m u d 




M n 


m u d 




M n 


0.02 




0.2927(18) 


0.03 




0.3570(26) 


0.04 




0.4081(21) 


m ud 


m s 


M K 


m ud 


m a 


M K 


m ud 


m 3 


M K 


0.02 


0.03 


0.3239(17) 


0.03 


0.02 


0.3288(28) 


0.04 


0.02 


0.3565(22) 


0.02 


0.04 


0.3528(17) 


0.03 


0.04 


0.3836(24) 


0.04 


0.03 


0.3830(20) 


0.02 


0.05 


0.3797(17) 


0.03 


0.05 


0.4088(24) 


0.04 


0.05 


0.4323(21) 



TABLE II: Scalar form factor /o(?max) a ^ simulated quark masses. 



m u d 


m s 


/(>(<?max) 


m u d 


m s 


/o(?max) 


m u d 


m a 


,/o(?Lx) 


0.02 


0.03 


1.00067(17) 


0.03 


0.02 


1.00050(22) 


0.04 


0.02 


1.00098(55) 


0.02 


0.04 


1.00202(48) 


0.03 


0.04 


1.00036(11) 


0.04 


0.03 


1.00024(10) 


0.02 


0.05 


1.00352(82) 


0.03 


0.05 


1.00126(35) 


0.04 


0.05 


1.00018(6) 



TABLE III: Fit parameters for interpolation of fo(q 2 ) to q 2 — using pole form Eq (13 1 fl . 





m 3 


X 2 /dof 


/o(0) 


A (1) 


0.02 


0.03 


0.76 


0.99955(47) 


1.16(39) 


0.02 


0.04 


0.67 


0.9979(14) 


1.13(32) 


0.02 


0.05 


0.60 


0.9952(25) 


1.09(28) 


0.03 


0.02 


1.49 


0.99872(43) 


2.23(40) 


0.03 


0.04 


1.91 


0.99912(26) 


1.75(32) 


0.03 


0.05 


2.07 


0.99697(90) 


1.60(30) 


0.04 


0.02 


1.98 


0.99574(95) 


1.97(33) 


0.04 


0.03 


1.79 


0.99911(20) 


1.80(28) 


0.04 


0.05 


1.06 


0.99924(15) 


1.61(23) 



14 



TABLE IV: Results for £(0) and / (0) obtained from Eqs. f32 ]l - 1[33 ]) . 









m s 




f(0) 






/o(0) 


0.02 
0.02 
0.02 






0.03 
0.04 
0.05 




-0.056(23) 
-0.101(34) 
-0.133(39) 






1.00007(94) 

0.9992(29) 

0.9962(57) 


0.03 
0.03 
0.03 






0.02 
0.04 
0.05 




+0.041(15) 
-0.0307(94) 
-0.053(16) 






1.00017(58) 
QQQ85('36 , 1 
0.9990(12) 


0.04 
0.04 
0.04 






0.02 
0.03 
0.05 




+0.050(22) 

+0.0195(85) 

-0.0144(65) 






0.9968(15) 

0.99915(31) 

0.99912(24) 


TABLE V: 


Fit parameters 


for chiral extrapolation of /+(0), and A/ at physical quark 


mass. 


fit form 




X 2 /dof 




CO 


Cl,s 


Cl,v 


A/ 


Eq. §7} 
Eq. (|38l 
Eq. 4MJ 




0.46 
0.46 
0.10 




-1.99(34) 
-1.2(2.0) 
-2.3(1.8) 


-44(14) 


-2.6(6.3) 
22.1(5.4) 


-0.013(2) 
-0.009(9) 
-0.003(11) 



TABLE VI: Recent lattice estimates of /+(0). Note that unquenched results in Ref.[34|,[35[ are preliminary. 
Two values in Ref. [34l] are obtained from two choices of the chiral extrapolation form (polynomial and ChPT 
based forms). 





Nf 


quark action 


/+(0) 


this work 


2 


domain-wall 


0.968(11) 


Becirevic et al. [33J 





improved Wilson 


0.960(9) 


JLQCD [34] 


2 


improved Wilson 


0.967(6), 0.952(6) 


MILC [35] 


3 


improved staggered 


0.962(11) 



TABLE VII: Fit parameter for chiral extrapolation of £(0). 



X 2 /dof 


di >v 


0.45 


-1.30(28) 
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FIG. 1: Effective mass plots for pion at m u[ j = 0.04 (top figure), 0.03 (middle figure), and 0.02 (bottom 
figure). Lines for data with zero meson momentum represent the fitted mass, while those for larger 
momentum are an estimation from the lattice dispersion relation Eq. (121[) . 



1G 



m J = 0.03, m = 0.05 




FIG. 2: Comparison of fitted energies (circles) with lattice dispersion relation Eq. (|21|l at (m u d,m s ) = 
(0.03, 0.05). Top and bottom panels show pion and kaon energies, respectively. 



m =0.03. m=0.04, t=10 




I i i i I i i i I I I I I I i i L_ 

20 40 60 80 

jack knife sample 



FIG. 3: Three-point functions (top and middle panels) and their double ratio R (bottom panel) for each 
jackknife sample at (m u d, m s ) = (0.03, 0.04). 
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FIG. 4: Double ratio Eq. ((22j) at m„ d = 0.04 (top figure), 0.03 (middle figure), and 0.02 (bottom figure). 
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FIG. 5: Plots of F(p, p') defined by Eq.pl]) with m s = 0.05 and |p'| = 0. Top and bottom figures show 
data with |p| 2 = l and 2, respectively. 
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FIG. 6: Plots of double ratio R k (t; p, p') defined by Eq. J26J with m s = 0.04 and |p'j =0. Top and bottom 
figures show data with |p| 2 = l and 2, respectively. 
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FIG. 7: Pole model interpolation of fo(q 2 ) to q 2 = 0. Top, middle and bottom figures show results at sea 
quark mass m u d — 0.04, 0.03 and 0.02, respectively. Filled circles are fo(q 2 ) at simulated q 2 , while open 
squares represent interpolated value to q 2 = 0. 
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FIG. 8: Interpolation of F(p, p') as a function of q 2 at m U( j = 0.03. Top and bottom figures show data 
with |p'| = and |p|=0, respectively. 
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FIG. 10: Chiral extrapolation of -Ra/- Top and bottom figures show results using Eqs, (|38[l and (1391) . 
respectively. 
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FIG. 11: Chiral extrapolation of £(g 2 ). 



